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Abstract. The implementation of a resistive- wall extension to the non-linear MHD-code 
JOREK via a coupling to the vacuum-field code STARWALL is presented along with first 
applications and benchmark results. Also, non-linear saturation in the presence of a resistive 
wall is demonstrated. After completion of the ongoing verification process, this code extension 
will allow to perform non-linear simulations of MHD instabilities in the presence of three- 
dimensional resistive walls with holes for limited and X-point plasmas. 



1. Introduction 

Plasma instabilities growing on fast time-scales such as vertical displacement events, disruptions, 
external kink modes, or edge localized modes are connected with a time-dependent magnetic 
perturbation outside the plasma. The mirror currents induced in conducting structures around 
the plasma 1 by the varying field perturbations act back onto the instabilities. Thus, linear and 
non-linear dynamics of the plasma can be strongly influenced by the wall. An external kink 
mode driven by strong edge current densities, for instance, is converted into a resistive wall 
mode (RWM) in the presence of a conducting wall close to the plasma edge and grows on the 
timescale of resistive wall-current decay as described in great detail in Refs. [1, 2]. 

We describe the status of a resistive wall extension for the non-linear MHD code JOREK [3, 4] 
as well as first benchmark results. The implementation is done by coupling JOREK to a modified 
version of the STARWALL code [5] which determines the magnetic field in the vacuum region 
in the presence of three-dimensional walls with holes. The "response" of conducting structures 
to magnetic field perturbations is represented by response-matrices computed by STARWALL. 
These matrices are used in JOREK to keep track of wall-currents and calculate the tangential 
magnetic field required for the modified boundary conditions as discussed in Section 3.2. 

The article provides some information about the STARWALL code in Section 2. The coupling 
equations are derived in Section 3. First benchmarks of the resistive wall extension are presented 

For simplicity, all conducting structures are denoted wall in the following. 



in Section 4. Conclusions and a brief outlook are given in Section 5. 



2. The STARWALL Code 

STARWALL solves the vacuum magnetic field equation outside the JOREK computational 
domain in presence of a three-dimensional conducting wall with holes as a Neumann-like 
problem. The continuity of the magnetic field component normal to the boundary of the JOREK 
computational domain (called interface) is used as a boundary condition. The wall is represented 
by infinitely thin triangles which allows to approximate realistic wall structures very well using 
an effective resistance rj w /d w where r? w and d w denote the resistivity of the wall material 2 and 
its thickness, respectively. Wall currents are assumed constant within each triangle such that 
they can be described by current potentials at the triangle nodes [5] 3 . 

In case of an ideally conducting wall, solving the vacuum field equation results in an algebraic 
expression for the magnetic field component tangential to the interface in terms of the normal 
component. Wall currents don't need to be considered explicitly as they are instantaneously 
given by the normal field at the interface. Using the ideal wall response-matrix M ld calculated 
by STARWALL, this can be written in the following way 5 : 
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Here, bi denotes a JOREK basis function consisting of a ID Bezier basis function along the 
poloidal direction on the interface and Fourier basis function in toroidal direction. denotes 
poloidal flux coefficients. The poloidal flux at the interface is given by $ = Ylj bj ^j- 

When a resistive wall is considered, wall currents cannot be eliminated anymore. In this case, 
the tangential magnetic field is given by 

= J> I^M- % + Y k J , (2) 
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where wall currents evolve in time according to 

n &&y h (3) 
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Here, M ee , M e y, M)" 5 , and denote resistive response matrices determined by STARWALL. 
Indices % and j run over all boundary degrees of freedom of the respective variable and k over 
all wall current potentials. For cross-checking, a relation between ideal and resistive response 
matrices can be derived by letting rj w — > oo (Appendix A). 

Information about the STARWALL code coupled with CASTOR (sometimes called 
STARWALL_C) can be found in Refs. [5, 7]. It is similar to the STARWALL code coupled with 
JOREK (sometimes called STARWALLJ). An article describing more details is in preparation 
by the STARWALL author Peter Merkel. 

3. Implementation in JOREK 

For the coupling of JOREK and STARWALL, Eqs. (2) and (3) are discretized in time as described 
in Section 3.1. The vacuum response then enters into a natural boundary condition in JOREK 
as discussed in Section 3.2. 

2 Wall resistivity is normalized the same way as plasma resistivity (see Ref. [6] for JOREK normalizations). 

3 For simplicity, we will (inaccurately) speak about wall currents when wall current potentials are referred to. 

4 Stan is the component of (B — B^) tangential to the interface, where B^ denotes the toroidal field component. 
For practical reasons, the algebraic expression is written in terms of the poloidal flux ^ instead of the normal 

field component (which can be calculated from of course). 



3. 1 . Time-Discretization 

Equation (2) is evaluated at the new timestep n + 1 and discretizations \I' n + 1 = ^ n + 5^> n and 
Y n+1 = Y n + 5Y n are used, where superscripts indicate evaluation at the given timestep. The 
tangential magnetic field is thus given by 
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The general time-evolution scheme described in Ref. [8], 
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is used for the JOREK equations written in the form <9X(u)/<9i = Z(u), where u denotes the 
vector of JOREK physical variables like temperature or density, X and Z the left- respectively 
right-hand side terms of the time-evolution equations, and 6 and £ are numerical parameters 6 . 
The same general time-evolution method needs to be applied to the wall-current evolution 
(Eq. (3)) giving 
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Terms with 5Y k are brought to the left hand side. The equation then reads 
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After solving for SYP, one gets 
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where some of the "derived response matrices" 

% bk, = mil 

A k>j = -(1 + Miys k , k Eij = M* + J2 k M% A k>j 
B Kk = -At^MH/S kik F hk = Mj(l + B k)k ) 
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6 Parameter values satisfying 8 — £ — 0.5 are required to guarantee second-order accuracy of the time-evolution 
scheme. For instance, (8 = 0.5, £ = 0) corresponds to a linearized Crank-Nicholson scheme and (8 — 1,£ = 0.5) 
to the linearized Gears scheme. 



have been used. These matrices are computed at the beginning of a JOREK simulation from the 
STARWALL response matrices. They need to be updated if one or more parameters entering 
the definitions have changed, e.g., At, £, or rj w . Plugging Eq. (8) into Eq. (4) gives 
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and, making use of Eq. (9) , we get 
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Boundary Integral 

The fixed boundary conditions for poloidal flux and plasma current corresponding to an ideally 
conducting wall are removed from JOREK. As a consequence, a boundary integral in the current 
definition equation resulting from partial integration (which vanishes in fixed-boundary JOREK) 
needs to be considered now. It can be written in terms of the tangential magnetic field such 
that we can insert Eq. (11) as a natural boundary condition. Details are given in the following. 
The current definition equation, j = A*ty , is written in weak form as 

/ dV §2 (■?' " A **) = °- ( 12 ) 

where j? denotes the test-functions taken to be identical with the basis- functions. Using 
A*^ = R 2 V • (R- 2 W), we get 

Jdvjtj-JdV j! V • (± = 0, (13) 

where the second term of the resulting expression can be integrated by parts ( j dV a V • b = 
— f dV Va • b + <f dA a b • n) yielding 

I dv W 2j + I dv W 2 VjT ' w " f dA R ( w ^/ R l = °- ( 14 ) 

= B t an 

Here, h denotes the unit vector normal to the interface and the tangential field is identified in 
the boundary integral (refer to Appendix B for details). Eq. (11) is inserted into Eq. (14) and 
after separating implicit and explicit terms, we have derived the form of the current equation 
implemented in JOREK, 
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Figure 1 

Left: JOREK flux-surfaces of a limiter equilibrium (limiter point shown as blue 
cross) are compared to CEDRES++ last closed flux-surface and magnetic axis. 
Right: Comparison of the respective q-prohles. 



with 
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where the integral over u and v is carried out by Gauss-quadrature: J du X(u) = 
Y^m G w m G X(u mc ) with weights w mc and positions u mG . Coordinate u corresponds to the 
element-local coordinate of the 2D Bezier elements oriented poloidally along the interface. The 
other element-local coordinate is denoted v. J2 denotes the 2D-Jacobian in the poloidal plane. 
Wall-currents are updated after each time-step according to Eq. (8) to guarantee consistency 
with the implicit time-stepping of JOREK. 

For the n = component, poloidal field coils need to be taken into account properly. This is 
done by replacing the expression for the tangential magnetic field, Eq. (2), by 

Stan = E 6 M E ^ (*i " *i° n ) + E Mil Y * + I ■ 
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Here, and -B£an i denote time-independent coil contributions to poloidal flux and tangential 
field at the interface, respectively. This leads to two additional terms in Eq. (15). A similar 
boundary integral also occurs in the Grad-Shafranov equation giving the plasma equilibrium. 
Here, only the n = component needs to be considered (axisymmetric equilibria in JOREK) 
and no time-discretization is required. 



4. First Benchmarks and Results 

In this Section, we present first results of the ongoing benchmark effort for code validation. 
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Figure 2 

Left: Linear growth-rates of a tearing-mode in a circular plasma with concentric 
ideal wall are compared between CASTOR using its own vacuum field module, 
CASTOR coupled with STARWALL, and JOREK coupled with STARWALL. Very 
good agreement is observed for a variety of different plasma resistivities and plasma- 
wall distances. 

Right: The current potential distribution (arbitrary units) reflecting the 2/1 tearing 
mode structure is plotted during the linear phase (rjj ore k = 1 • 10~ 6 , r w = 1.2). 



4-1. Free-boundary equilibrium 

Computing a free-boundary equilibrium with JOREK and comparing it to the results of 
another equilibrium solver allows to check some parts specific for n = like poloidal field- 
coil contributions (see Sec. 3.2). We are considering an ITER- like limiter plasma. As seen 
from Figure 1, flux surfaces and q-profile of the JOREK equilibrium agree very well with the 
equilibrium computed by the CEDRES++ code [9]. Slightly different discretizations of the 
poloidal field coils might be responsible for small remaining differences. 

4-2. Tearing Mode 

A 2/1 tearing mode in a circular plasma with major radius R = 10, minor radius a = 1 and 
uniform plasma resistivity surrounded by an ideally conducting wall is considered. Figure 2 
shows excellent agreement between JOREK and the linear CASTOR code [10]. Repeating 
the simulation for a resistive wall with zero resistivity is physically identical, of course, but 
numerically treated in a different way. The exact agreement we get here proves consistency. 

4.3. Resistive Wall Mode 

In this Section, results for 2/1 resistive wall modes (RWMs) in a circular plasma with minor 
radius a = 1 and major radius R = 10 surrounded by a resistive wall are presented. The left 
part of Figure 3 shows linear growth rates obtained for different wall radii r w as a function of the 
wall resistivity. Clearly, the modes cannot be stabilized beyond a certain wall distance from the 
plasma even at very low wall resistivities. The right part of the Figure shows energy time-traces 
during the linear and saturation phase of a typical simulation. 
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Figure 3 

Left: Linear RWM growth rates extracted from the simulations are shown for a 
variety of wall radii and wall resistivities. 

Right: Energy time-traces are shown for linear and saturation phase of an RWM 
(r w — 1.2, n wa u = 1 • lO -4 ). In saturation, wall currents on the high-held side have an 
amplitude comparable to the low-held side while they are smaller in the linear phase. 



5. Conclusions and Outlook 

The ongoing implementation and verification of a resistive wall model in the non-linear MHD- 
code JOREK was summarized. Benchmarks for a freeboundary equilibrium and tearing mode 
cases show good agreement. Simulations of the linear and non- linear phase of RWMs were 
presented (comparison to analytical theory and linear codes ongoing). 

Benchmarking will be continued and extended to nondinear comparisons. Realistic X-point 
geometries will be considered requiring a special treatment of grid corners. After completion, 
the code can be applied to a variety of MHD instabilities interacting with conducting structures 
like resistive wall modes, edge localized modes, vertical displacement events or disruptions. 
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Appendix A. Relation between ideal and resistive response matrices 

Between the ideal and resistive response-matrices, the relation 



k 

must hold which becomes obvious when letting r/ w — > and inserting Eq. (3) into Eq. (2). Also, 
the no-wall response can easily be identified as 

MfJ = M% (A.2) 

when letting Yj. — > 0. Eqs. (A.l) and (A.2) can be used to cross-check some parts of the 
implementation as both forms of the wall-treatment with eliminated wall-currents (Eq. (1)) and 
explicitly treated wall-currents (Eqs. (2) and (3)) are implemented in JOREK. 



Appendix B. Magnetic field tangential to the interface 

The magnetic field in the JOREK reduced-MHD model is given by 
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B = — e^ + - W x e^. 
K K 
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Its component tangential to the boundary of the JOREK computational domain (in the poloidal 
plane) can be determined from 

Btan = (Bxn)-e^ (B.2) 

where n denotes the unit normal vector of the interface and e^ the normalized toroidal basis 
vector. Inserting Eq. (B.l) and using vector identities, this can be written in the form 
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which can be identified in the boundary integral of Eq. (14) and replaced by the STARWALL 
response. 



